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The diagnosis of wind energy conversion systems (WECS) turns out to be 
necessary because of their relatively high cost of operation and maintenance. 
Wind turbines are hard-to-access structures, and they are often located in 
remote areas. Therefore, a remote diagnosis (e-diagnosis) is required. This 
paper proposes an alternative approach for the e-diagnosis of a WECS based 


on the discrete wavelet transform (DWT) and frequency analysis of the aero 
generator stator currents. To validate this approach, real-time hardware in the 
Keywords: loop (HIL) is used to simulate in real-time the mathematical model of the 
induction generator on the OPAL-RT OP5600 platform to generate the stator 
currents and the rotor speed. The DWT is applied to the current signal, to 
. generate the DWT signal, which has a huge number of points that are not 
Hardware in the loop supported for direct transmission by the Arduino Mega RobotDyn because of 
Wavelet transforms its limited sample time. The absolute values of the DWT peak points (MDWT) 
Wind turbine are sent as point’s packages form to the diagnosis station via the ESP8266 
integrated Wi-Fi board of the Arduino Mega RobotDyn to monitor the SCIG 
states and determine the number of broken bars. 
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1. INTRODUCTION 

Wind energy conversion systems (WECS) are becoming increasingly important throughout the world, 
due to the exponential decreasing of fossil fuels and increasing environmental concerns such as global warming 
effect and pollution [1], [2]. Nowadays, many wind turbines and industrial fields use asynchronous squirrel 
cage motors. Since they are robust, simple to build, easy to use and require minimum low-cost maintenance. 
However, as robust as they are, they keep being exposed to many failures that conduct, after a certain time, 
partial damage or a complete stop. The presence of a fault can trigger emergency shutdowns in an industrial 
process. Which causes losses in production processes [3]. To avoid these undesirable situations, fault detection 
and diagnosis methods, that help to detect faults since their first appearance, were invented [4]—[6]. 

The importance of early detection of fault locations and time in electric generators has caught the 
researchers’ attention because it prevents the spread of damage and increases generator life [7]—[11]. Wavelet 
transform (WT) is a technique of signal analysis based on the use of region windows of variable size. Because 
the wavelets have a limited duration and frequency bandwidth, the WT is localized in the frequency as well as 
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time domains. The signal is represented by some coefficients that make the WT a powerful tool for both fault 
diagnosis and condition monitoring [12], [13]. 

The automated control of modern wind machines requires proactive maintenance [14] because of the 
operating and maintenance costs and due to the remote location of wind turbines and their difficult structures, 
remote diagnosis (e-diagnosis) is necessary. Research in this field is still new and only a few have been 
mentioned, one of which is a method where data is collected by various sensors installed on each wind turbine 
and transmitted via wireless communication to the monitoring station [9]. Subramaniyan and Ramesh [10], the 
authors used GPRS/GSM communication and embedded system integration with ZigBee wireless sensor 
network technology for monitoring and remote control. 

This work aims to develop a centralized architecture based on the use of the OPAL-RT OP5600 
platform to simulate in real-time the mathematical model of the squirrel cage induction generator (SCIG) and 
generate stator currents and rotor speed signals. The discrete wavelet transform (DWT) is applied to the current 
signal and the abstract values of the DWT output signal peak points are sent out to the Arduino Mega RobotDyn 
to be transmitted via Wi-Fi to the diagnosis station to monitor the SCIG states. 

Among the problems of the e-diagnosis, encountered in this experiment, is the limitation of the 
Arduino board because of the high noise in the analog input measurements and the Wi-Fi transmission 
frequency. Thus, in this newly developed strategy, the DWT is applied on packages of 20ms input data to 
reduce the number of points transmitted at a lower frequency while maintaining the precision of the information 
as will be detailed later. The proposed e-diagnosis strategy has been experimentally validated in real-time 
hardware in the loop (HIL) simulation and has shown good results in fault monitoring and surveillance. 

The article is organized as follows: in the 2" section, the model of WECS is described; the principle 
of DWT application is presented in the 3“ section; in the 4" section, the proposed e-diagnosis strategy of the 
WECS is presented and validated in the experimental simulator; the experimental results and discussions are 
presented in the 5" section, and the conclusion is drawn in the 6" section. 


2. MODELING OF THE WIND ENERGY CONVERSION SYSTEM 
The basic diagram of the WECS is presented in Figure 1. The wind turbine drives a SCIG equipped 
with a capacitor excitation system. The models of these elements are presented as follows. 


Gear Box 


Fixed capacitor bank 


i 
Blades Squirrel Cage 
1 


Induction Generator 


Figure 1. SCIG with a capacitor excitation system driven by a wind turbine 


2.1. Model of the wind turbine 
Wind turbine blades are used to convert wind kinetic Energy into electrical energy via the coupled 
generator. The mechanical power captured by the wind turbine is expressed by (1) [15], [16]. 


Pa = = pAv3 C(A, p) (1) 


Where, p is the air density in Kg/m’, A is the area swept by the turbine blades in m? and is equal to 7R?, where 
R is the radius of the wind turbine blade, v is the wind speed in m% and C,(A, p) is the power coefficient. The 
C, (A, P) is a non-linear function of the tip speed ratio A and the pitch angle f, which is given by (2) and (3) [17]. 
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aes wrR (4) 
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Where, wr is the rotational speed of the blades (rads). 


2.2. The SCIG stator equations 
For more than three decades, many mathematical models-based techniques have been created [18]- 
[21]. The stator is given by (5). 


; dQs 
Vs = Rsis + Te (5) 
Were 
Vs = [VaVpVel" (6) 
is = liaipie]" (7) 
7,00 
R; = [ozo (8) 
00r, 
Øs = Leis + Lsrir (9) 


Lms Lms 
Ls = -3 Lis + Ls — ~ 10) 


Lis and Lms are leakage and magnetizing inductances of the stator windings, respectively. 


Lms = #2 N? (3) a1) 


The matrix of mutual inductance stator-rotor is given by (12). 


LarLaz «Lan, 
Loy = LyiLp2 . Lon, (12) 
Lerkce Leon, 


Where: 
LNs |. r r 

Loi = ar sin 5 cos (9; +“) (13) 
Or: the i” bar angle, a,: the angle between two adjacent bars. 
2.3. Model of the healthy rotor 

Figure 2 presents the equivalent circuit of the rotor cage, containing k+1 magnetically meshes. The 
rotor's bars are numbered from 1 to k+1. The two neighboring bars are joined by a segment end ring to define 
each mesh. In addition, as shown in Figure 2, each end ring segment and rotor bar is represented by a resistor 


and inductor equivalent circuit. The rotor equation can be written: 


y dór 
Vr = Rii, + AE (14) 
With: 


9, = Liris + Lt, (15) 
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j TAA : ka Pe 
i, = iste RE eo ae tendo (17) 
T 
0, = [9,19,2 we Ørk e Ørn, Del (18) 
The rotor resistance matrix ( N, x N, ): 
R, —R, O .. 0 -Rp 
—R, Rog —R, =. 0 
RSE eas AET 
. ; Ro — Rp 
—=Rp . . m — Rp Ro (19) 
Ro =2(Rp + Re) (20) 
Lo = Lmr + 2(Lp + Le) (21) 
lr r 
Line = "2" a, (1-7) (22) 
Lo Iy2-Lp l a Lyw,-1) Lin, — Ly 
Lo — Lp Lo L23 — L, r. . Lon, 
=| ba bea Ly bo 
Ly, = Lp . . wae LN,(Np-1) = Ly Lo (23) 
_ _ Holr (a 
La = — "2" (2) (24) 
The expression of the torque is given by (25). 
pt OLsr(Or) ; 
D = 5 Pig i, (25) 


P: rotor poles pair number, 0,: the electric displacement angle of the rotor in radians 
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Figure 2. Equivalent circuits of a SCIG healthy rotor 


2.4. Rotor modeling with rotor fault 

As shown in Figure 3, the number of rotor equations is reduced according to the broken bars number, 
and the meshes damaged are eliminated since the imesh will be n times larger. The mutual inductance matrix 
stator rotor with dimension 3x(N, —n) becomes: 


Lar Laz ERD Lai Laçi+(n+1)) KEN Lan, 
Ler = |Lo Loz =- Lpi Loser) - Low, (26) 
Lea Lez ios Lei Leçirin+) Dos Len, 
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After eliminating n rows and columns of rotor matrix inductance, the new matrix of (N, — n) x 
(N; — n) dimension becomes: 


Lo Iyg-—Lyp Lii Ly (i4n41) B Lin, — Lp 
L21 — Lp Lo = Lzi Lz(i+n+1) aa Lon, 
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Ly,1 — Lp Ly,2 ise Ly, i Ly, (i¢n41) ee Lo (28) 
Loi = Lii + 2(Lp + (n + 1)Le) (29) 
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The rotor resistance matrix with (N, —n)x(N, —n) dimension becomes: 
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Figure 3. Equivalent circuits of a SCIG rotor with n broken bars 


3. DISCRETE WAVELET TRANSFORM 

The SCIG current signal faults, which appear as oscillations in the current signal spectrum, are 
obtained using the current signature analysis approach. The choice of the mother wavelet is completely 
arbitrary. As a result, selecting the mother wavelet correctly reduces errors in fault diagnosis [22], [23]. 

The temporal current signals are represented in the time-frequency spectrum using the symlets-DWT 
signal processing approach to detect the WECS SCIG broken bars number. Therefore, by selecting an initial 
scale to acquire the current signal characteristics, several representations of the time-frequency spectrum with 
different features can be observed. The augmentation of this scale allows to obtain more detailed information 
in the time-frequency spectrum. 

The family of wavelets is obtained from a mother wavelet in Figure 4 denoted Ys (t) (t = s = 1) 
by expansion-translation (by varying t and s). This mother wavelet is in the form: 


1 t-s 
bsd =F) (34) 
With s and 7 are its scale and translation respectively and s, t E IR and T # 0 [24], [25]. 
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Figure 4. Symlets mother wavelet signal 


4. PROPOSED FAULTS E-DIAGNOSIS STRATEGY OF WECS 
The proposed strategy of WECS fault e-Diagnosis, as illustrated in Figure 5, goes through the 
following three main modules: 


4.1. RT-LAB | OP5600-OPAL-RT 
In this module, five operations are executed: 

— HIL simulation of the SCIG and wind turbine assembly to generate stator currents and rotor angular speed 
signals (the used sample time in the HIL simulation is 0.3 ms). 

— Application of the symlet DWT, order 8, levels 1 on packages of 64 points of the current signal to provide 
a processing signal technique for data packages of approximately 20 ms to examine and detect, in real-time, 
any perturbations and quick frequency variations in the current signal. 

— After correlation between the current signal, obtained from de HIL simulation, and the mother wavelet, the 
current signal will be decomposed at different levels called detail and approximation coefficients. A 
comparison between the healthy currents and the failed ones (one, three, and five broken bars of SCIG) 
with different fault levels is performed by the display of coefficients which are transformed into the scale 
translation coefficients. To be able to take a decision on the fault’s diagnosis, a simplification of the function 
DWT is proposed as (35): 


MDWT (t,) = |DWT(t,)| (35) 
With txis the instances where eee 
The MDWT output signal is transmitted as an electrical signal between 0 and 4.5 V, via the FPGA analog 
output port of the OP5600-OPAL-RT. 


= 0, suchas 0 < t < tpwith tp ~ 20ms is the package time range. 


4.2. Arduino RobotDyn 
In this case, three operations are executed: 
— For the e-diagnosis, an Arduino RobotDyn control board with a built-in ESP8266EX Wi-Fi module is used 
to send real-time data from RT-Lab to the diagnosis station in a remote wireless way. 
— The electrical signal generated by RT-Lab is captured by the analog input of the Arduino and sent in 8 bits 
with a speed of 115200 baud to the Wi-Fi module. 
— The Wi-Fi module starts sending data once it is connected to an access point. 


4.3. Diagnosis station 

A PC equipped with MATLAB/Simulink is connected to the same wireless access point as the 
Arduino and starts receiving data via TCP/IP protocol. A MATLAB/Simulink program is developed for real- 
time data receiving and showing the generator’s health. The program implemented in the diagnosis station can 
identify in real-time the number of broken bars based on the MDWT variations. 
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Figure 5. Proposed diagnosis flowchart 


5. HIL EXPERIMENTAL RESULTS AND DISCUSSION 


Oo 


1091 


For the real-time HIL simulation and e-diagnosis of the isolated WECS, the experimental prototype 
shown in Figure 6 is used while Table 1 shows the main prototype components. As demonstrated in Figure 7, 
the WECS model and the current DWT are implemented in the OPAL-RT OP-5600 RT-Lab platform that 
generates the instantaneous MDWT signal via the analog output port. This signal is transmitted in real-time to 
the e-diagnosis station via Wi-Fi using an Arduino Mega data acquisition and transmission with a built-in 
Wi-Fi board where technical characteristics are presented in Table 2, while characteristics of the diagnosis 
station are presented in Table 3. The characteristics of the access point that interconnects the e-diagnosis station 
and Arduino are given in Table 4. The SCIG parameters used in the HIL real-time simulation for the different 


operating conditions are given in Table 5. 


Table 1. Components of the developed experimental HIL simulation prototype 


N° Designation 


OP5600 real-time digital simulator 


TP-LINK Wi-Fi N 300Mbps Access Point 
Host PC equipped with RT-LAB software 
Power Bank 

Oscilloscope 

e-diagnosis station 


ADUNABRWNKE 


Arduino compatible Mega R3 ATmega2560 and Wi-Fi ESP8266 with 32MB memory 
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Figure 6. Photo of the experimental test bench Figure 7. Block diagram of the proposed e-diagnosis strategy 


Table 2. Technical characteristics of the RobotDyn Arduino 
Parameters Description 
Microcontroller ATmega2560 (16 Mhz)+ESP8266 (Esp8266 connects to rx3/tx3 of atmega2560) 
Flash memory 32Mb, 


Interface Micro-USB-TTL CH340G(converter ch340g USB connects to rx0/tx0 of atmega2560) 

Pin input/output Digital I/O Pins 70, 16 analog inputs, 14 PWM, 4 UART buses. Output current for 5 V - about 800 mA, for 
3.3 V — 180 mA 

Pigtail connector Connecting an antenna on this device, you will increase your reach, right from 90 meters to 240 meters 

Supply voltage Operating voltage 5 V 


Table 3. Diagnostic station characteristics 


Parameters Description 
Processor Intel(R) Core (TM) i5-3337U CPU @ 1.80 GHz 
Memory Installed 8192 MB RAM 
System type 64-bit operating system, Windows 10 Pro. 
MATLAB version 2018b 
Hard disk SSD (solid state drive) 500 GB 
Wireless Wi-Fi 802.11 b/g/n 


Table 4. Access point characteristics 


Parameters Description 
Interface 10/100M Ethernet port (RJ45) 
Wi-Fi standards IEEE 802.11n, IEEE 802.11g, IEEE 802.11b 
Frequency 2.4-2.4835 GHz 
Signal rate 11n: Up to 300 Mbps (dynamic), 11g: Up to 54 Mbps (dynamic), 11b: Up to 11 Mbps (dynamic) 
Transmit power <20 dBm (EIRP) 
DHCP DHCP Server 


Table 5. WECS parameters 


Parameters Symbol Values Parameters Symbol Values 
Rated power (kW) P 4 Rotor bar resistance (uQ) Rp 96.94 
Rated voltage (V) V 220/380 Rotor bar inductance (uH) Ly 0.28 
Rated current (A) I 15.2/8.8 | End-ring segment resistance (uQ) Re 5 
Rated speed (rpm) Q 1465 End-ring segment inductance (uH) Le 0.036 
Pole-pair number p 2 Number of rotor bars Nr 28 
Supply frequency (Hz) fs 50 Mean radius air-gap (mm) r 54 
Power factor coso 0.83 Effective air-gap (mm) g 0.28 
Rotor inertia (kg.m°) J 0.011 Length of air-gap (mm) l 120 
Friction coefficient (Nm/rad/s) f 8.73e-4 Wind(m/s) v 7 
Stator phase turns in series N, 156 Radius of the wind turbine blade(m) R 1.2 
Stator phase resistance (Q) Rs 1:5 Blades pitch angle(deg) B 8 
Stator phase leakage inductance (mH) Lls 7 Air density (Kg/m*) p 1.225 
Stator magnetizing inductance (H) Lms 0.55 
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To verify the efficiencies of the generalized model, the control strategy, and the diagnosis procedure, 
three simulation cases of the isolated WECS shown in Figure 1 are validated considering that the SCIG of the 
WECS is subjected to one, three, or five broken bars. In the real-time simulation, the defects are introduced by 
removing the blocks that define the broken bars and inserting new inductances and resistances as indicated in 
section 2. The suggested approach was tested on a4 kW WECS with the detailed parameters given in Table 5. 
The WECS real-time simulation using HIL has been implemented with a sample time of T, = 0.3 ms. 

From Figure 8 or Figure 9, it is difficult to notice the change in current variations for one or three broken 
bars respectively, but when five bars are broken, ripples can be observed in the stator current signal shown in 
Figure 10. For one broken bar speed signal in Figure 11, it is hard to figure out the speed changes, compared with 
3 and 5 broken bars in Figure 12 and Figure 13 respectively where the speed decreasing can be observed. 


pas a E Eerie pm E a 


= M 100s 
M Pos:640,ams M Pos:640.0ms 
Figure 8. Stator current for 1 broken bar Figure 9. Stator current for 3 broken bars 


ES r _ eee ES ewe 


M 160s TEE ; 
M Pos:640,0ms M Posi-270s 
Figure 10. Stator current for 5 broken bars Figure 11. Rotor angular speed for 1 broken bar 


i) E pm ee 


5.005 = ESA = 
M Pos:648.8ms M Pos:640,8ms 


Figure 12. Rotor angular speed for 3 broken bars Figure. 13. Rotor angular speed for 5 broken bars 


However as shown in Figure 14, Figure 15, and Figure 16, which present the DWT signal, as well as 
Figure 17, Figure 18, and Figure 19, which present the MDWT signal, it is easy to notice the variations of 
DWT and MDWT in the instances of bars breakings. Figure 20 shows the real-time MDWT signal received by 
the diagnosis station for 1, 3, and 5 broken bars. The signal values extracted from this figure allow us to know 
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the number of the broken bars, while Table 6 indicates the reference interval of MDWT values to determine 
the number of broken bars. This can be used to find an indirect relation between the MDWT values and the 
number of broken bars in the SCIG rotor. 


M1005 a 
M Pos :6.00ps. M Pos :8.08ps 


Figure 14. DWT signal for 1 broken bar Figure 15. DWT signal for 3 broken bars 


ME re r rere 


100s CHi= 2,00) M500s 
M Pos:6.08ps M Pos:640.0ms 
Figure 16. DWT signal for 5 broken bars Figure 17. MDWT for 1 broken bar 


SS e SS eee 


penne 


CHi= 1,60) 
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Figure 18. MDWT signal for 3 broken bars Figure 19. MDWT signal for 5 broken bars 


5 broken bars 


3 broken bars 


20 40 60 80 100 120 140 160 180 200 220 
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Figure 20. MDWT variation for 1, 3 and 5 broken bars 


Int J Pow Elec & Dri Syst, Vol. 14, No. 2, June 2023: 1085-1097 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 O 1095 


Table 6. MDWT values for different rotor broken bars 


MDWT interval Number of broken bars 
0 < MDWT < 6- 107° 0 
6- 1073 < MDWT < 0.5 1 
0.5 < MDWT < 1.5 3 
MDWT > 1.5 5 


6. CONCLUSION 

The paper proposes a novel strategy for an isolated wind energy conversion system (WECS) 
e-diagnosis based on the discrete wavelet transform (DWT) and frequency analysis of the aero generator stator 
currents. This e-diagnosis strategy has been implemented in real-time for an asynchronous squirrel cage 
induction generator (SCIG) based on the DWT (symlet, order 8, levels 1) of the stator current signal to 
determine the generator's broken bars number. This proposed strategy is validated by a real-time Hardware in 
the loop (HIL) simulation of the mathematical model of the WECS SCIG using the RT-Lab OPAL-RT OP5600 
platform. The DWT is applied to the generated current signal, it has a huge number of points which makes its 
transmission by Arduino Mega RobotDyn impossible because of the limited sample time of this later. The 
absolute values of the DWT peak points (MDWT) are sent as points packages via the Arduino integrated 
ESP8266 Wi-Fi module to the diagnosis station to monitor WECS states and determine the number of SCIG 
broken bars based on the MDWT variations. The obtained results demonstrate the applicability of this 
technique for the real-time analysis of different electrical variations and the instant detection of different faults. 
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